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Abstract 



We systematically describe and classify 1-dimensional Schrodinger equations that can be solved 
in terms of hypergeometric type functions. Beside the well-known families, we explicitly describe 2 
new classes of exactly solvable Schrodinger equations that can be reduced to the Hermite equation. 

1 Introduction 

Exactly solvable 1-dimensional Schrodinger equations play an important role in quantum physics. The 
best known are the harmonic oscillator and the radial equation for the hydrogen atom, which are covered 
in every course of quantum mechanics. A number of other examples were discovered in the 30's of the 
last century. They include (known also under other names) the Eckart potential [Eck , the trigonometric 
and hyperbolic Poschl- Teller potencjal |PT| and the Rosen potential [MRj . All these examples can be 
reduced to the hypergeometric equation (see eg. |WW| ) . One should also mention the Morse potential 
that leads to the confluent equation. Problems involving these potential are often used in classes on 
quantum mechanics, see eg. the well-known problem book of Fliigge [Fiji] . A number of techniques have 
been developed to study their properties, such as the factorization method HI] and the method of the 
superpotential (see for instance |CKSj ). 

In later years J. N. Ginocchio [Gin] discovered that these examples can be generalized to a larger 
class of potentials equivalent to the hypergeometric or confluent equation. Later this class was extended 
by G. A. Natanzon [Nat] and further generalized by R. Milson [Mil] . These classes, besides the hitherto 
known potentials, are not very practical in applications, since they are not given by explicit expressions. 

There exist even more general classes of potentials that can be called exactly solvable, found by 
A. Khare and U. Sukhatme |KSj . They are, however, expressed in exotic special functions. 

In our paper we will systematically describe basic classes of exactly solvable Schrodinger equations. 
Most of the material of our paper is scattered in the literature, notably in }Natj and }Mil] (see also |Bos] ) . 
Our treatment is somewhat more systematic than what one can find in the literature. For instance, we 
distinguish between the complex classification and the real classification. It also seems that for the first 
time we explicitly describe 2 new classes of exactly solvable Schrodinger equations, which can be reduced 
to the Hermite equation (see Subsects 15.21 and I5.3[) . 

By a (stationary 1-dimensional) Schrodinger equation we will mean an equation of the form 



An operator of the form (|1.2p will be called a (1-dimensional) Schrodinger operator. V{r) will be called 
a potential and the parameter E an energy. 

(|1.2|) can be interpreted as an operator in a number of ways. If V is a holomorphic function of the 
complex variable on some open VL C C, then r can be interpreted as a complex variable and (|1.2p can be 
viewed as an operator on holomorphic functions. In this case it is natural to allow E to be a complex 
parameter. The corresponding eigenvalue equation (|1.1[) will be then called a Schrodinger equation of 
complex variable. Let us note that, if a ^ 0, b are complex constants, a complex affine transformation 
r i->- ar + b preserves the class of Schrodinger equations of the complex variable. 




(1.1) 



This equation can be interpreted as the eigenvalue problem for the operator 



H := -d 2 r +V{r). 



(1.2) 
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One can also interpret r as a real variable in some open Jet The operator (II. 2[) is then viewed as 
acting on functions on /. Understood in this way, will be called a Schrodinger equation of the real 

variable. Of special interest is then the case of real V . 

Our paper is organized as follows. First, we will briefly dicuss some general facts related to 2nd 
order linear differential equations. In particular, we will describe basic ingredients of the so-called _Bo.se- 
Natanzon method, which permits to obtain a class of Schrodinger equations equivalent to a chosen linear 
equation. It is worthful to point out, that it applies to more general situations than those described in 
the literature ( [Nat] . [Mil] ) . We will formulate a criterion which determines when the Bose-Natanzon 
method can be used. 

Next, we will focus on hypergeometric type equations. Recall that equations of the form 

(a(z)d 2 + r(z)d z+ r 1 )f(z)^0, (1.3) 
where o~(z), t(z) are polynomials with 

degcr < 2, degr < 1, 

and i] is a number, are occasionally called hypergeometric type equations (see A.F. Nikiforov and V.B. Uvarov 
|NU) ) . Solutions of (|1.3[) are called hypergeometric type functions. Hypergeometric type functions are 
very well understood and include the hypergeomeric function, the confluent function, Hermite, Laguerre 
and Jacobi polynomials, etc. Traditionally, one classifies hypergeometric type equations into several dis- 
tinct classes invariant under complex affine transformations of the form z > az + b (o ^ 0). In each class 
one chooses a simple representative, to which the whole class can be reduced. Such representatives are 
the hypergeometric equation, the confluent equation, etc. 

In our analysis, we will concentrate on 3 specific equations: 

(1) the hypergeometric equation 

(z(l - z)d 2 + (c-(a + b+ l)z)d z - ab) f(z) = 0; (1.4) 

(2) the generalized confluent equation 

(zd 2 z + (c- 1 z)d z -a) f(z) = 0, (1.5) 

(essentially just the rescaled confluent equation); 

(3) the generalized harmonic oscillator 

(-d 2 z + 6 2 z 2 + pz + A) f{z) = 0, (1.6) 

(essentially just the translated harmonic oscillator equation). 

Note that all hypergeometric type equations with a nonzero a are related with one of the above 
equations by elementary operations (multiplication of f(z) and of both sides of the equation by a function) 
and by a change of variables not depending on the parameters a, 6, c (respectively a, c, 7 or l9, /?, A). 

(1) If a is 2nd order and has 2 distinct roots, an affine transformation reduces (|1.3|) to (|1.4|) . 

(2) If a is 2nd order and has only one root, then a transformation involving z } z^ 1 reduces (|1.3[) to 
(TOl) with 7^0. 

(3) If a is 1st order, then an affine transformation reduces (| 1 . 3[) to (jl.5p . 

(4) If <7 is Oth order and nonzero, then (|1.3p is equivalent to (|1.6p with (9^0. 
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Obviously, (|1.4I) and (II .5|) are themselves equations of the hypergeometric type. (|1.6p is not, but can 
be reduced to a hypergeometric type equation. This reduction is depends on whether 6 ^ or 9 = 0. 
The former case leads to a hypergeometric type equation with a constant er, as mentioned above. In the 
case 8 — 0, (|1.6|) is the Airy equation, which can be reduced to a special case of the hypergeometric type 
equation by a transformation involving z^z 3 . 

Equations (|1.4j) . (|1.5j) and (jl.6j) are very well understood and have well known solutions. We will 
describe classes of solvable potentials reducible to one of the equations (|1.4|) . (|1.5j) or (|1.6j) . We will 
describe both the complex and the real classification of such potentials. 

All solvable potentials that we consider are holomorphic on the complex plane, apart from some 
isolated singularities. Therefore, it is natural to consider first a classification of Schrodinger equations of 
the complex variable. Obviously, we can always use a complex affine transformation to put the equation 
in a convenient form. One can also move a complex constant from the potential to the energy. 

Suppose we are given a Schrodinger equations of the complex variable. Embedding a real line (or its 
subset) in the complex plane in an arbitrary way yields a Schrodinger equation of the real variable. In 
order to obtain an interesting and useful real classification we will demand in addition that the potential 
is real. Thus, we will look for an open, possibly infinite interval / C R, on which the potential is real 
for a selected class of parameters. It turns out, that in all cases that we consider J, can be chosen so 
that it ends either at — oo, +oo or at a singularity of the potential, and no singularities lie inside /. Such 
an interval / will be called a natural real domain. Again, a real affine transformation and moving a real 
constant from the potential to the energy can be used to put the equation in a convenient form. 

Clearly, if J is a natural real domain, then the operator H = + V(r) is hermitian on C£°(I). It is 
then natural to ask about self-adjoint extensions of H . We will not discuss this question here. We plan 
to consider it in later papers. 

In Equations (|1.4I) . (| 1 . 5|) and (jl.6|) we have 3 arbitrary parameters. Therefore, in all the solvable cases 
that we will describe, the potential depends on 2 parameters, since the third parameter is responsible for 
the energy E. 

Acknowledgement The research of J. D. was supported in part by the grant N N201 270135 of the 
Polish Ministry of Science and Higher Education. 

2 Second order homogeneous differential equations 

The main object of our paper are ordinary homogeneous 2nd order linear differential equations, that is 
equations of the form 



We will then say that the equation (|2.7p is given by the operator (|2.8p . 

We will treat r either as a complex or a real variable. In the complex case we will usually assume 
that the coefficients are analytic. 

In this section we describe some general facts related to the theory of equations of the form (|2.7j) and 
their reduction to the Schrodinger equation. 

2.1 The Bose invariant of a second order differential equation 

Suppose that we are given a second order differential equation (12.7[) . Then we can always eliminate the 
first order term as follows. We divide from the left by a(r), set 




(2.7) 



It will be convenient to treat (|2 . T[) as the problem of finding the kernel of the operator 

A(r,d r ) := a{r)d 2 r + b{r)d r + c(r). 



(2.8) 
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We check that 

h{r)A{r,d T )h{r)- x = d 2 +I(r), 
where, using the notation a' := -P-a(r), 



iac - 2ab' + 2bo! - b 2 
Alp 



l =— ■ ( 2 - 9 ) 



Clearly, <p{r) := h(r)<f>(r) solves the equation 

(dl+I{r))4>{r) = 0. (2.10) 

(|2.10p will be called the canonical form of the equation (12.71) . (|2 .9[) will be called the Bose invariant 
of (|2.7p [Bos]. (The name "Bose invariant" is used eg. by Milson [Mil] . however the object itself was 
clearly known before Bose). 

It is easy to check that for any functions /, g the equation 

f{r)A(r,d r )g(r)4,(r)=0 (2.11) 



has the same Bose invariant as (|2.7p . Conversely, if an equation has the same Bose invariant as (|2.7[) , it 
equals (|2.1ip for some /, g. 

2.2 Schwarz derivative 

The Schwarz derivative of r(y) appears naturally in the context of 2nd order ordinary differential 
equations and was known to XIX century mathematicians. It was mentioned already by Kummer in 1836 
[Kumj . The univaersally adopted name Schwarz derivative was introduced by A. Cayley, who derived its 
basic properties in 1880 [Cay]. For more historical comments together with numerous applications we 
refer the reader to a book by Osgood |Osg| . 
The Schwarz derivative is defined as 



r'{y)J 2\r'(y) 
Here are some basic properties of the Schwarz derivative: 
Proposition 2.1 1. Let y H ► s(y), s r(s) be two functions. Then 

{r,y} = {r,s} S '(y) 2 + { S ,y}. (2.13) 

2. Let r t— > y(r) be the inverse of y i— > r(y). Then 

{y,r} = -{r>(y)]- 2 {r,y} (2.14) 

3. {r,y} = 0^r(y) = ^,ad-bc^0 
We will need the following fact: 

Proposition 2.2 Suppose that z' = z' p (\ — z) q . Then 

{z, r} = ^z 2p - 2 {l - z) 2q - 2 ((p 2 - 2p)(l - z) 2 + {q 2 - 2q)z 2 - 2pqz{\ - z)) . (2.15) 
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2.3 The Liouville transformation 

Consider a 2nd order equation in the canonical form, that is 

(<9 r 2 + /(r))0(r) =0. (2.16) 
Let us make the transformation r — r(y) in this equation. We obtain 

- T d 2 r cj>(r(y)) ~ T7T^drHr(y)) + /(r(tf))0(r(y)) = (2.17) 



[r'(y)] 2 r " v yy " [r'(y)] 3 
The resulting equation we transform again to its canonical form: 

(d 2 + J(y))i>(y) = 0, (2.18) 

where 

^y) = 7^]2<^(y))' ( 2 - 19 ) 

J(y)^[r , (y)] 2 I(r(y)) + ^{r,y}, (2.20) 

The above procedure is called the Liouville transformation of \2. 1 6}) by the change of variables r = r(y) 
[Lioj . One can check using (I2.17[) , that the composition of two Liouville transformations, first by the 
change of variables r = r(y), then y = y(z), is the Liouville transformation by the change of variables 
r = r{y{z)). 



2.4 Natanzon's problem 

Let A(di, ■ ■ . , a.k', z, d z ) be a 2nd order differential operator depending on k parameters, such that its 
Bose invariant can be written as 

I(b 1 ,...,b n ;z) = b 1 I 1 (z) + ... + b n I n {z) (2.21) 

for some (at least one) n < k linearly independent functions I\, . . . , I n (independent of a±, . . . , a^), and 
some n numbers bi(a\, . . . , a^). Then one can solve the following problem (which we will call Natanzon's 
problem for A(z,d z )): 

Find all potentials V(r) such that the 1- dimensional stationary Schrodinger equation 

{-d 2 r + V(r) - E) <P{E, r) = (2.22) 

can be transformed to the equation given by A(a%, . . . , a^; z, d z ) for some a\(E), . . . , ak{E). We allow the 
following operations: 

(1) multiplication of both sides of the equation by some f(E,r); 

(2) substitution of (f)(E,r) := g(E,r)(j)(E,r) for some g(E,r); 

(3) change of coordinates r i— > z independent of E. 



This is the problem solved by Natanzon in the case of the hypergeometric equation and by Milson in the 
general case of hypergeometric type equations. In the following, we recall the construction they used, in 
a slightly generalized form. 

Let us consider an arbitrary equation of the form ()2.22|) . Obviously, its Bose invariant equals E—V(r). 
Clearly, the transformations (1) and (2) allow us to transform A(a\, . . . , a^; z, d z ) to its canonical form 

-d 2 z -I(bx,...,b n ;z). (2.23) 



G 



Thus (I2.22[) can be transformed to A{a-y, . . . , a/-; z, d z ) if 

I(h, . . . , b n ; z) = (r'(z)f (E - V (r(z)) + l{r, z}, (2.24) 

that is, if the two Bose invariants are related by a Liouville transformation. Using (|2.2ip we rewrite 
(|2~23|l as 

b x {E)h{z) + ... + b n {E)I n {z) = {r'{z)f (E-V(r{z)) + ^{r,z}. (2.25) 

By assumption, 6, depend on E and h(z) are independent of E. Thus, the dependence of 6, on E is 
linear. Therefore, by transforming linearly (b\, . . . , b n ) into (&i, . . . , 6 n _i, E), we can assume that 

I(bi, . . .,b n ;z) = 61/1(2) + . . . + 6 n _i/ n _i(z) + EI n (z), 

for some functions Ii(z), which are linear combinations of Ii(z). Then we can write 

EI n (z) + b 1 I 1 (z) + ... + b n _ 1 I n _ 1 (z) = (r'(z)) 2 (E-V(r(z)) + ^{r,z}. (2.26) 
Therefore, Natanzon's problem is solved by the following pair of equations 

V(h, 6 n _i,r(z)) := (r'(z)y 2 (bj^z) + ... + 6 B _ 1 7„_ 1 (z) - i{r, z} 

(r'(z)f = I n {z). (2.27) 
We can rewrite (|2.27l) in an equivalent way as 

V(pi, . . . , 6„_i,r) := (z'(r)) 2 (6 1 / 1 (z(r)) + ... + &„_ 1 /„_ 1 (z(r))) +i{z,r}, 

(z'(r)y 2 = I n (z(r)). (2.28) 

If desired, we can renormalize E by substracting a constant. 

Note that in practical applications the value of the above method depends on whether it is possible 
to invert the relation r <— > z expressing it in terms of standard functions. In the next sections, we 
will describe such solutions of the Natanzon problem for the three classes of equations described in the 
introduction. 

3 Schrodinger operators reducible to the hypergeometric equa- 
tion 

3.1 Hypergeometric equation 

The hypergeometric equation is given by the operator 

F(a, b; c; z, d z ) := z(l - z)d\ + (c-(a + b + l)z)d z - ab, (3.29) 

where a, b, c are arbitrary complex parameters. Its Bose invariant equals 

. , , (1 - a 2 - b 2 - 2ab)z 2 + 2(-2ab + ac + be - c)z + 2c- c 2 
2h(a,b;c;z) = 4 z 2 (l - z) 2 ' 
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(4)77277— ( 5 )l72> (6)777— 72- (3-30) 



It satisfies condition (|2.21[) for n = 3. The functions Ii(z), i — 1,2,3 can be taken to be 4z ^_ z ) 

STiTTP ano - 4z(i-z) • ^ * s na tu r al to demand that (r'(z)) 2 = ^(z) is a function of the form z 2p (l — z) 2q . 
It is easy to see that this gives the following possibilities: 

(1) 4z(l-z)' (2) 4z 2 (l-z)' (3) 4z(l-z) 2 ' 

4: 2 {l^ W (5) 4^' (6) 4(T^z 

We will see that each of these ansatzes leads to exactly solvable potentials considered in the literature. 
Note also that the formula for the Schwarz derivative (j2.15[) will be handy. 

3.2 The Riemann equation 

It is well known that it is useful to consider the hypergeometric equation as a special case of the so-called 
Riemann equation. The Riemann equation is defined as the class of equations on the Riemann sphere C 
having 3 regular singular points. The following theorem summarizes the basic theory of these equations: 

Theorem 3.1 (1) Suppose that we are given a 2nd order differential equation on the Riemann sphere 
having 3 singular points Zi,z 2 ,Z3, all of them regular singular points with the following indices 

zi ■ Pi, Pi, 

Z2 ■ P2,P~2, 

Z3- P3, Pi- 
Then the following condition is satisfied: 

Pi +pl +P2 +P2 + P3 +P3 = 1- (3.31) 

Such an equation, normalized to have coefficient 1 at the 2nd derivative, is always equal to 

Zl z 2 z 3 

V I Pi P2 P3 z,d z I <j>{z) =0, (3.32) 



Pi P2 P3 



where 



zi z 2 z 3 



vipi P2 P3 z, dz : = a 2 2 -(^^ + ^ r i + ^ i 



Pi P2 P3 

I PlPl(zi-Z 2 )(z 1 -Z 3 ) , P2P2(Z2-Z 3 ){Z 2 -Z 1 ) , p 3 p 3 (z 3 - Z 1 ) (z 3 - Z 2 ) 
(Z — Zl) 2 (z-Z2)(z-Z3) (Z-Z2) 2 (Z-Z3)(Z — Zl) (z - Z 3 ) 2 (z - Z 1 ) (z — Z2 ) ' 

(2) Let z i— > w = h(z) = g^xg ■ (Transformations of this form are called homographies or Mobius 
transformations). We can always assume that ad — be — 1. Then 

h(zi) h{z 2 ) h(z 3 ) \ / zi z 2 z 3 

V \ Pi P2 P3 w,d w ^(cz + d^vi pi P2 P3 z,d z 



(3) 



Pi P2 P3 I \ Pi P2 P3 



Zl Z 2 



( z - zj) A (z- z 2 y 



Zl 


Z2 Z 3 


Pi 


P2 P3 


Pi 


P~2 P3 


Z3 


z,d z \ 


P3 




P3 
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Clearly, in all above formulas one of can equal oo, with an obvious meaning of various expressions. 
For convenience we give the expression for the Riemann operator with 2:3 = 00: 

Z\ z 2 00 
V I Pi Pi P3 z,d z 
Pi h P3 

= #-( Pl + *- 1 + P2 + i *- v \e M 

(3.33) 



z — Z\ z — z 2 

PlP~\{zx - Z 2 ) p1pl{Z2 - Zl) p-ip-i 



(z - zi) 2 (z - z 2 ) {z - z 2 ) 2 (z - zi) (z - zi)(z - Z2) 

The hypergeometric equation is a special case of the Riemann equation, since 

/ 1 00 

J-(a,b;c; z,d z ) = z(l — z)V \ a z,d z 

I 1-c c — a — b b 

= z(l ~ z)dl + (c~{a + b+ l)z)d z - ab. 

Every Riemann equation can be brought to the form of the hypergeometric equation by applying (1) and 
(2) of Theorem O 

From Theorem 13.11 we also see that symmetries of the hypergeometric equation are better visible if 
we replace a, 6, c with a, /?, fi: 

a = c — 1 (3 = a + b — c, fj, = a — b; 
a= i+g±g-fl 6= I±££±/3±H c=l + a. 

In fact, the new parameters coincide with the differences of the indices of the points 0, 1, 00: 

a = Pi- Pi, P = P2 - P2, P = P3 - Pa- 
in the parameters a, /3, /i, the Bose invariant of the hypergeometric equation has a more symmetric form: 

(l-q 2 )(l-z) + (l-/? 2 )z + ( M 2 -l)(l-z)z 
2h{a,p,Mz) = 4z2(1 _ ^2 ■ 

We can summarize the relation between the hypergeometric equation and its canonical form by 



- z)^- x T l + a; Z) fl,) - z)- 



2 2 



-V 



3.3 Trigonometric Poschl- Teller potentials 

In this subsection we consider Ansatz (1) of (|3.30p . We set 

z = sin 2 r, which solves z' = 2z? (1 — z)i. (3.35) 



9 



This leads to the Schrodinger equation 

(-d? + V a ,p(r) ~ V 2 ) Hr) =0, (3.36) 

where 

4 / sin r V 4 / cos r 

?2 1 \ /( /„,2 o2\ 



/3 2 1\ 4 (a 2 -/? 2 )4cos2r 



4/ sin 2 2r 2 sin 2 2r 



This potential was proposed and solved by G. Poschl and E. Teller [PTj . It is usually called the Poschl- 
Teller potential, sometimes also the Poschl- Teller potential of the first kind. 

A natural real domain for this potential is ]0, 2tt[. The potential is real if r €]0, 27r[ and a 2 ,/? 2 are 
real. 

Explicitly, the reduction of (|3.36|) to the hypergeometric equation is derived as follows: 



-4^+3(i - z )f +3 t (s+e+B+i, z+i-E+i. i + a - Zi dz ) 2-* -i (i - zy 

1 oo 

-4z(l-z)V\ f + i f + i f z,d z 

_a , 1 i 1 _£i 

2 4 2 ' 4 2 



2 4 



= -4 z (l-z)\d 2 z +(±-- 1 



2z 2(1 -z) 
4 7 z ' V 4/1 



^ 2 -])-+(^-i)t^-^ (3-38) 



^+(- 2 -l)- i 2 -+(/ ?2 -l)^ 2 - 
V 4/sin 2 r V 4/cos 2 r 



3.4 Hyperbolic Poschl- Teller potentials 

We continue with (1) of ([33D]) . We set 

z=-sinh 2 r, which solves z' = -2(-z)i (1 - z)5 . (3.39) 
This leads to the Schrodinger equation 

(-^ + V^(r)-/i 2 )^>) =0, (3.40) 



where 



2 M 1 L.3 1 



4 / sinh 2 r \ 4 / cosh 2 r ^ ^ 

2 +/3 2 1\ 4 (a 2 -^ 2 )4cosh2r 



4/ sinh 2 2r 2 sinh 2 2r 



This potential was also proposed and solved by G. Poschl and E. Teller }PTj . In the literature it is known 
as hyperbolic, generalized Poschl- Teller potential or the Poschl- Teller potential of the second kind. 

A natural real domain for this potential is ]0, oo[. The potential is real if r (E]0, oo[ and a 2 , 1 are real. 

To see that (|3.42p can be solved in terms of the hypergeometric equation, we first repeat the compu- 
tations leading to (|3.38l) . and then set z = — sinh 2 r. 

Note that the trigonometric and hyperbolic Poschl- Teller potentials are related with a complex linear 
transformation r > ir. Therefore, these two potentials are two real cases of the same complex case. 
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3.5 Alternative approach to Poschl- Teller potentials 



Ansatzes (2) and (3) lead to the same classes of exactly solvable potentials as Ansatz (1). To see this it 
is enough to consider Ansatz (3). In fact, we first repeat computations analogous to Q3.38[) : 



1 oo 

-4^(1 - zfV I f + i | f + i z,d z 



2 1 4 2 2 1 4 

Oil _§_ _{± I 1 
2^ 4 2 2 T 4 



^-iJ^+^rb-K-j.i <'-■■■' 



We set 



2 * 

2 = tgh r, which solves z = 2z 3 (l — z), (3.42) 
obtaining the trigonometric Poschl- Teller potentials, or 

z = -tg 2 r, which solves z' = -2(-z)' (1 - z), (3.43) 

obtaining the hyperbolic Poschl- Teller potentials. 

3.6 Rosen/Scarf/Woods-Saxon potentials 

Let us consider Ansatz (3) of (|3.30p . We set 

z = y~ 9 r -. which solves z' = 2z(z—l). (3.44) 

We also set n = (3 2 — a 2 . This leads to the Schrodinger equation 

{-d 2 r + K, M (r) + a 2 ) 0(r) = 0, (3.45) 

where 



V Kttl (y) := n- 



M 2 1 



4lKy! ' 2coshr V 4 47 cosh 2 r 

This potential was proposed and solved by M. F. Manning and N. Rosen |MRj . In the literature it is also 
called the Scarf [Sea) or Woods-Saxon [SWj potential (for instance in |Fhi] ) . 

A natural real domain for this potential is ] — oo, oof. The potential is real if r G] — oo, oo[ and k, j3 2 
are real. 

Here is a derivation of (|3.45p from the hypergeometric equation: 

-4z 1+ t(l - z) 1 +#^( °+^+ 1 , i + a . z ,d x )z-*(l - z)-$ 

= -Az 2 (z-l) 2 V 
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= -4z 2 (l-z) 2 ^ + +a 2 (l-z)+/3 2 z-( M 2 -l)z(z-l) (3.46) 

= -a 2 + a 2 + (- a 2 +/3 2 ) T ^ : -( / , 2 -l)^^. (3.47) 

3.7 Eckart/Hulthen potentials 

We still consider Ansatz (3) of (|3.30p . 

We set 

z = - which solves z' = 2z(l — z). (3.48) 

We also set n = f3 2 — a 2 . This leads to the Schrodinger equation 

(~d 2 + V K ^(r) + a 2 ) tf>(r) = 0, (3.49) 

where 

e r ffi 2 1\ 1 



V K Av) : - K 2sinhr ' V 4 4 J sin h 2 



r 



This potential was proposed and solved by C. Eckart [Eckj . In the literature (for instance [Flu ) it is 
also called the Hulthen potential [Hulj . sometimes also the generalized Morse potential, because of its 
similarity to the Morse potential, see Subsect. 14.41 

A natural real domain for this potential is ]0, oo[. If r G]0, oo[ and k, (3 2 are real, then the potential is 
real. 

To check these computations, we first repeat the computations leading to (|3.46l) . and then set z = 
l 



Note that from the Rosen and Eckart potentials are related with a complex affine transformation 
r H> m — r. Therefore, these two potentials can be treated as two real cases of the same complex case. 

3.8 Alternative approach to Rosen/Scarf/Wood-Saxon and Eckart/Hulthen 
potentials 

Ansatzes (5) and (6) lead to the same classes of exactly solvable potentials as Ansatz (4). To see this it 
is enough to consider Ansatz (5). In fact, we first repeat computations analogous to (13. 46[) : 

-4z-^(l - z) 1+ ^( a+0+ / +1 , 9±S±JH1 ] i + a - Zl d z )z-^(\ - z)~$ 




-4(i - zf (a 2 - j^- z d 2 

+ {a 2 - l)^—^- + P 2 - - [i 2 — -. 

z* z z 



We set 



= l-e 2r , which solves z' = 2(z - 1), (3.50) 
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obtaining the Rosen potential, or 

z = l + e 2r , which solves z' = -2(z - 1), (3.51) 
obtaining the Eckart potential. 

4 Schrodinger operators reducible to the confluent equation 

One of basic exactly solvable equations is the confluent equation, given by the operator 

T(a; c; z, d z ) := zd 2 z + (c - z)d z - a. (4.52) 

It is convenient to consider (|4.52l) in parallel with the equation given by the operator 

F(a, b;-;z, d z ) := z 2 d 2 z + (-1 + (1 + a + b)z)d z + oh. (4.53) 

The equation given by (I4.53|) is sometimes called the 2-F0 equation. Note that 

z a T(a, b; — ; z, d z )z~ a = wJ-(a; 1 + a — b; w, d w ), w = —z , z = —w . 

Hence the 2-F0 equation is equivalent to the confluent equation. The relationship between the parameters 
is 

c = 1 + a — 6, 6 = 1 + a — c. 

Another exactly solvable equation that we will consider in this section is sometimes called the 0-F1 
equation. It is given by 

T(c;z,d z ) := zd 2 + (c~ z)d z -a, (4.54) 

and is equivalent to the Bessel equation. 

Clearly, the confluent, 2^0 and qFi equations belong to the class of hypergeometric type equations. 
The basic equation of this section will be (|1.5j) . given by 

F(a; c; 7; z; d z ) := zd 2 z + (c — r yz)d z - a. (4.55) 

It will be called the generalized confluent equation. Note that in the case 7=1 (|4.55j) coincides with the 
confluent equation. If 7 7^ 0, (|4.55[) can be reduced to the confluent equation by rescaling (and hence 
also to the 2-F0 equation). If 7 = 0, (|4.55|) coincides with the 0-F1 equation. 
The Bose invariant of the generalized confluent equation equals 

T , . , -7 2 z 2 + 2(c7-2a)z + 2c-c 2 
1/1 (a, b; 7; z) = 



Az 2 

4.1 Symmetries of the generalized confluent equation 

Let us first describe symmetries of the generalized confluent equation. 

z c - 1 T{a;c;"/;z 1 d z )z 1 ~ c = T(a - c~f + 7; 2 - c; 7; z, d z ); (4.56) 
e^ 7Z J r (a; c; 7; z, d z )e JZ = — T{c — a; c; 7; w, d w ), z — —w; (4-57) 

Besides, the scaling acts as follows: 

F(a; c; 7; z, d z ) = 7^(0/7; c; 1; w, d w ), w = 72. 
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It is convenient to introduce new parameters a, v: 

a = ' 2 ' — , c = 1 + a: 

a = c — X = a — 6, v — cy — 2a = 1 — a — b. 

In the new parameters the Bose invariant of the generalized confluent equation has a more symmetric 
form: 

j 2 v 1 — a 2 

Thus the starting point for the further analysis will be the equation 

l! _ JL + (— - -) — 

4 2z \ 4 4/z 2 



^ + -T-sr+ T-7 3 UW = 0. (4.58) 



We will describe 3 classes of Schrodinger operators solved using the confluent equation corresponding 
to three obvious choices for (r'{z)**j : 

(1) \, (2) ± (3) -L. (4.59) 

4.2 Hydrogen atom 

We consider (1) of (|4.59|) . We set 

z = 2r. (4.60) 

This leads to the Schrodinger equation 

{-d 2 r + V a , u (r) + 7 2 ) 0(r) = 0, (4.61) 

with the potential 

T , v fa 2 l\ 1 

(|4.61|) is the radial part of the Schrodinger equation for the hydrogen atom. 

]0, oo[ is a natural real domain for this equation. If v, a 2 are real, then so is V aM on ]0, oof. 
The derivation of (|4.6ip from (|4.58[) is immediate: 

4(-d 2 + ^-- + (^--)- 
V 4 2z V 4 4J z 2 

4.3 Rotationally symmetric harmonic oscillator 

We consider (2) of ([435]) . We set 

z = r 2 , which solves z' = Isfz. (4.62) 
This leads to the Schrodinger equation 

{-d 2 r +V a , 1 {r)-2v)4>{r) = 0, (4.63) 
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with the potential 



]0,oo[ is a natural real domain. For real j 2 ,a 2 , the potential is real on ]0,oo[. By scaling, the only 
diffcret real cases are 7 2 = 1,0,-1. 

For 7 2 = 1 the equation is the radial part of the rotationally symmetric harmonic oscillator. For 
7 2 = the equation is the radial part of the Helmholtz equation. 

Here is an explicit derivation of (|4.63[) from (|4.58|) : 



i / , 7 2 v (a 2 1\ 1 

4/ z 



z l 



= -Azd 2 z - 2d z + 7 2 z -2v + [a 2 - 
= -9 2 +7 V-2, + ( a 2 -i)l. 



4.4 Morse potential 

We consider (3) of 1)435]) . We set 

— T 

z = c , 

This leads to the Schrodinger equation 



which solves z 1 = —z. (4.64) 



a 2 



-d 2 + V^(r) + —) ( f>(r)=0, (4.65) 
where 

K, 7 (r) = ^e- 2 ''-^-'-. 

] — oo,oo[ is a natural real domain. For real v, 7 2 , the potential is real on ] — oo, oo[. By translation, 
the only different real cases are 7 2 = 1, 0, —1. 

Here is an explicit derivation of (|4.65l) from (|4.58l) : 

^-M-* + ^-i + f^-l)^i«« 

z 2 d 2 — zd z 



2 

7 Z v 


+ 


(-■ 


T ~ Yz 




V 4 




V 


a 2 


2 


+ T 




+ 


a 2 


2 




4 ' 



5 Schrodinger equations reducible to the generalized harmonic 
oscillator 

The last family of exactly solvable Schrodinger equations that we describe is (|1.6[) . given by 

- d 2 + 9 2 z 2 + pz + A. (5.66) 
We will call (|5.66|) the generalized harmonic oscillator equation. It is already in the canonical form. 
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If 8 2 7^ 0, then (|5.66l) is the Schrodinger equation for translated harmonic oscillator. It is then 
equivalent to the Hermite equation, given by 



G(a,y,d y ) :=(d 2 -2yd y ~2a). 



In fact, we have 



8e " 2 G(a,y,d y )e v 2 
8 + y 2 + 2a - 1) 



,2 



where y = s/9 (z + . 

For = 0, p 7^ 0, (|5.66p is equivalent to the Airy equation 

(d 2 + y)1>(y) = 0, 

which in turn is equivalent to a special case of the qFi equation: 

—d 2 + pz + A 

= P H-dl + y ) 
21 (2 

= (p3)3 W 3j-( -; W ,dy 

where we set y = pi z + p^^\, and w = 3~ 2 y 3 . 

We will describe 3 classes of Schrodinger operators that can be reduced to (|5.66p corresponding to 



three obvious choices for 



(r'(z))' 



(1) 1, (2) z, (3) z 2 . (5.67) 



5.1 Generalized harmonic oscillator 

We consider Ansatz (1) of (|5.67l) . which corresponds to the most obvious choice of the energy that is —A. 
We rename the variable z = r. This leads to the Schrodinger equation 

{-d 2 + V e , p (r) + A) <j>{r) - 0, (5.68) 

where 

V e , p (r) =9 2 r 2 +pr. 

5.2 Special potential I 

We can choose the energy to be —p, that is Ansatz (2) of (|5.67p . This corresponds to the substitution 

)2 
3 j 
. which solves z' = z~?. 

and leads to the Schrodinger equation 

(-d 2 + V e .Ar)+p)Hr) = 0, (5.69) 
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where 



In fact, 



2 ' 



2 



zi- 1 + 9 2 z 2 + pz + X) z~i 

1 „9 1 „ 9 A / 1 1\ 1 

- -*- 2 (!) § — (D s -a4)(i) 

Note the following intriguing feature of the above potential: the coefficient at r~ 2 is fixed and one 
cannot change it by rescaling the variable r. 

5.3 Special potential II 

We can choose the energy to be — 9 2 , that is Ansatz (3) of (|5.67[) . This corresponds to the substitution 

z = (2r) 5 ; which solves z' = z -1 , 
and leads to the Schrodinger equation 

(-d 2 + V p ,x{r)+9 2 )^r) = 0, (5.70) 



where 



In fact, 



, , p A 3 1 
(2r)a 2r 16 r z 



z5- 2 {-d\ + 9 2 z 2 +pz + X) z~i 

2 ' 2 2 



±d 2 + ± dz+ 9 2 + P + ± ( l ' ^ 1 



1 1 



(2r)i 2r \2 2 2 J 2 2 r 2 ' 
Again, the coefficient at r -2 is fixed and one cannot change it by rescaling the variable r. 
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Abstract 

We systematically describe and classify 1-dimensional Schrodinger equations that can be solved 
in terms of hypergeometric type functions. Beside the well-known families, we explicitly describe 2 
new classes of exactly solvable Schrodinger equations that can be reduced to the Hermite equation. 



1 Introduction 

Exactly solvable 1-dimensional Schrodinger equations play an important role in quantum physics. The 
best known are the harmonic oscillator and the radial equation for the hydrogen atom, which are covered 
in every course of quantum mechanics. A number of other examples were discovered in the 30's of the 
last century. They include (known also under other names) trigonometric and hyperbolic Poschl- Teller 
potentials |PT] , the Scarf potential |Scaj , the Eckart potential |Eckj , the Manning- Rosen potential MR 
and the Rosen-Morse potential. All these examples can be reduced to the hypergeometric equation (see 
eg. [WW]). One should also mention the Morse potential that leads to the confluent equation. Problems 
involving these potentials are often used in classes on quantum mechanics, see eg. the well-known problem 
book of Fliigge [Flu] . A number of techniques have been developed to study their properties, such as the 
factorization method |HIj and the closely related method of the superpotential (see for instance [CKS 
and [Cot]). 

In later years J. N. Ginocchio [Gin] discovered that these examples can be generalized to a larger 
class of potentials equivalent to the hypergeometric or confluent equation. Later, this class was extended 
by G. A. Natanzon [Nat] and further generalized by R. Milson Mil] . These classes, besides the hitherto 
known potentials, are not very practical in applications, since they are not given by explicit expressions. 

There exist even more general classes of potentials that can be called exactly solvable, found by 
A. Khare and U. Sukhatme [KS| . They are, however, expressed in exotic special functions. 

The literature on exactly solvable Schrodinger equations is very large. The subject is in fact very 
useful for applications, especially in quantum physics. We also believe that it is quite beautiful, which in 
the existing literature is perhaps not so easy to see. 

In our paper we would like to systematically describe basic classes of exactly solvable Schrodinger 
equations. Most of the material of our paper is scattered in the literature, notably in [Nat] . [Mil] and 
[CKSj (see also (Bosj ). Our treatment is, however, somewhat more systematic than what we could find in 
the literature. For instance, we distinguish between the complex classification and the real classification. 
It also seems that for the first time we explicitly describe 2 new classes of exactly solvable Schrodinger 
equations, which can be reduced to the Hermite equation (see Subsects 15.21 and 

By a (stationary 1-dimensional) Schrodinger equation we will mean an equation of the form 

{-dl + V(r) - E) <t>{r) = 0. (1.1) 

This equation can be interpreted as the eigenvalue problem for the operator 

H := -df,+V(r). (1.2) 

An operator of the form (11.21) will be called a (1-dimensional) Schrodinger operator. V(r) will be called 
a potential and the parameter E an energy. 

(|1.2p can be interpreted as an operator in a number of ways. If V is a holomorphic function of the 
complex variable on some open C C, then r can be interpreted as a complex variable and (|1.2[) can be 
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viewed as an operator on holomorphic functions. In this case it is natural to allow E to be a complex 
parameter. The corresponding eigenvalue equation (jl.ip will be then called the Schrddinger equation of 
the complex variable. Let us note that, complex affine transformations preserve the class of Schrodinger 
equations of the complex variable. (By a complex affine transformations we mean r i— > ar + b, where 
a 0, b are complex constants). 

One can also interpret r as a real variable in some open I cl. The operator (|1.2j) is then viewed as 
acting on functions on I. Understood in this way, will be called the Schrodinger equation of the real 
variable. Of special interest is then the case of real potential. Clearly, real affine transformations preserve 
the class of real Schrodinger equations of the real variable. 

Our paper is organized as follows. First, we will briefly dicuss some general facts related to 2nd 
order linear differential equations. In particular, we will describe basic ingredients of the so-called Bose- 
Natanzon method, which permits to obtain a class of Schrodinger equations equivalent to a chosen linear 
equation. Note that it applies to more general situations than those described in the literature ([Nat], 
Mil ). We will formulate a criterion that determines when the Bose-Natanzon method can be used. 

Next, we will focus on hypergeometric type equations. Recall that equations of the form 

(o-(z)d 2 z +T(z)d z + rj)f(z) = 0, (1.3) 
where cr(z), t(z) are polynomials with 

degcr < 2, degr < 1, 

and 77 is a number, are sometimes called hypergeometric type equations (see A.F. Nikiforov and V.B. Uvarov 
|NUj ). Solutions of (| 1 .3[) are called hypergeometric type functions. Hypergeometric type functions are very 
well understood and include the hypergeomeric function, the confluent function, Hermite, Laguerre and 
Jacobi polynomials, etc. Traditionally, one classifies hypergeometric type equations into several distinct 
classes invariant under complex affine transformations. In each class one chooses a simple representative, 
to which the whole class can be reduced. Such representatives are the hypergeometric equation, the 
confluent equation, etc. 

In our analysis, we will concentrate on 3 specific equations: 

(1) the hypergeometric equation 

(z(l - z)d 2 z + (c-(a + b+ \)z)d z - ab) f(z) = 0; (1.4) 

(2) the rescaled confluent equation 

(zd 2 + (c- 7 z)d z -a)f(z) = 0. (1.5) 

(3) the translated harmonic oscillator 

(-d 2 + 6 2 z 2 + pz + \)f(z)=0. (1.6) 

Note that all hypergeometric type equations with a nonzero a are related with one of the above 
equations by elementary operations 

• multiplication of f(z) by a function, multiplication of the equation by a function, 

• a change of variables not depending on the parameters a, 6, c (respectively a, c, 7 or 6, p, A). 

(1) If a is 2nd order and has 2 distinct roots, an affine transformation reduces (| 1 . 3[) to (1 1 .4[) . 

(2) If a is 2nd order and has only one root, then a transformation involving z i-> z _1 reduces (|1.3[) to 
(fL5)) with 7^0. 
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(3) If a is 1st order, then an affine transformation reduces (|1.3[) to (|1.5[) . 

(4) If a is Oth order and nonzero, then (|1.3p is equivalent to (|1.6I) with # ^ 0. 

Obviously, (|1.4|) and ()1.5j) are themselves equations of the hypergeometric type. (|1.6j) is not, but can 
be reduced to a hypergeometric type equation. This reduction is depends on whether 8 ^ or 8 = 0. 
The former case leads to a hypergeometric type equation with a constant c, as mentioned above. In the 
case 8 = 0, (1 1 . 6|) is the Airy equation, which can be reduced to a special case of the hypergeometric type 
equation by a transformation involving z n- z 3 . 

Equations Q1.4p . (jl.5[> and (jl.6|) are very well understood and have well known solutions. We will 
describe classes of solvable potentials reducible to one of the equations (|1.4I) . (|1.5[) or (|1.6[> . We will 
describe both the complex and the real classification of such potentials. 

It is natural to consider first a classification of exactly solvable Schrodinger equations of the complex 
variable. Note that all solvable potentials that we consider are holomorphic on the complex plane, apart 
from some isolated singularities. Obviously, we can always use a complex affine transformation to put the 
equation in a convenient form. One can also move a complex constant from the potential to the energy. 

Consider a family of exactly solvable holomorphic potentials. Suppose that for a selected subfamily of 
parameters the potential is real if restricted to an open, possibly infinite interval JcC. By elementary 
properties of holomorphic functions, if we extend the interval /, the potential is still real until we hit 
a singularity. A real affine transformation can be used to put the equation in a convenient form. In 
particular, we can always assume that / is a subset of the real line. We can also move a real constant 
from the potential to the energy. 

The above discussion motivates the following definition: We will say that an open interval I C M. is 
a natural real domain for a certain family of potentials if it ends either at — oo, +oo or at a singularity 
of the potential, and no singularities lie inside /. When describing the real classification, we will always 
specify a natural real domain of the potential. 

Clearly, if / is a natural real domain, then the operator H = —df, + V(r) is hermitian on C£°(J). It is 
then natural to ask about self-adjoint extensions of H . We will not discuss this question here. We plan 
to consider it in later papers. 

In Equations (|1.4j) . (|1.5|) and (|1.6|) we have 3 arbitrary parameters. Therefore, in all the solvable cases 
that we will describe, the potential depends on 2 parameters, since the third parameter is responsible for 
the energy E. 

Acknowledgement The research of J. D. was supported in part by the grant N N201 270135 of the 
Polish Ministry of Science and Higher Education. 

2 Second order homogeneous differential equations 

The main object of our paper are ordinary homogeneous 2nd order linear differential equations, that is 
equations of the form 

(a(r)d 2 + b(r)d r + c(r)) <j>{r) = 0. (2.7) 
It will be convenient to treat (|2.7[) as the problem of finding the kernel of the operator 

A(r,d r ) := a(r)d 2 r + b(r)d r + c(r). (2.8) 

We will then say that the equation (|2.7p is given by the operator (|2.8[) . 

We will treat r either as a complex or a real variable. In the complex case we will usually assume 
that the coefficients are analytic. 

In this section we describe some general facts related to the theory of equations of the form (|2 . T[) and 
their reduction to the Schrodinger equation. 
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2.1 The Bose invariant of a second order differential equation 

Suppose that we are given a second order differential equation (|2.7p . Then we can always eliminate the 
first order term as follows. We divide from the left by a(r), set 

Mr):=eXP (/ ^T) dt 
We check that 

where, using the notation a' :— ^a(r), 

4ac - 2ab' + 2ba' - b 2 , . 

/= 4^ ■ (2 - 9) 

Clearly, cf>(r) :— h(r)4>(r) solves the equation 

(d 2 + I(r)) for) = 0. (2.10) 

(|2.10[) will be called the canonical form of the equation (|2.7p . (12.91) will be called the Bose invariant 
of (|2.7|) [Bos]. (The name "Bose invariant" is used eg. by Milson [Mil] . however the object itself was 
clearly known before Bose). 

It is easy to check that for any functions /, g the equation 

f{r)A(r,d r )g(r)4,(r)=0 (2.11) 



has the same Bose invariant as (12.71) . Conversely, if an equation has the same Bose invariant as (12.7[) , it 
equals (|2.11[) for some /, g. 

2.2 Schwarz derivative 

The Schwarz derivative of r h-> r(y) appears naturally in the context of 2nd order ordinary differential 
equations and was known to XIX century mathematicians. It was mentioned already by Kummer in 1836 
[Kum] . The universally adopted name Schwarz derivative was introduced by A. Cayley, who derived its 
basic properties in 1880 [Cay]. For more historical comments together with numerous applications we 
refer the reader to a book by Osgood [Osg| . 
The Schwarz derivative is defined as 

{T y r= mx-i(m\ ( , 12) 



r'(y)J 2\r'(y) 
Here are some basic properties of the Schwarz derivative: 
Proposition 2.1 1. Let y i— > s(y), s i— > r(s) be two functions. Then 

{r iy } = {r, S } S '(y) 2 + { Sl y}. (2.13) 

2. Let r i— > y(r) be the inverse of y i— > r(y). Then 

{y,r} = ~[r'(y)r 2 {r,y}. (2.14) 

3. { r , y }^0^r(y) = ^,ad-bc^0 
We will need the following fact: 

Proposition 2.2 Suppose that z' — z p (l — z) q . Then 

{z, r} = ^z 2p - 2 {l - z) 2q - 2 ((p 2 - 2p)(l - z) 2 + {q 2 - 2q)z 2 - 2pqz{\ - z)) . (2.15) 



5 



2.3 The Liouville transformation 

Consider a 2nd order equation in the canonical form, that is 

(<9 r 2 + /(r))0(r) =0. (2.16) 
Let us make the transformation r — r(y) in this equation. We obtain 

- T d 2 r cj>(r(y)) ~ T7T^drHr(y)) + /(r(tf))0(r(y)) = (2.17) 



[r'(y)] 2 r " v yy " [r'(y)] 3 
The resulting equation we transform again to its canonical form: 

(d 2 + J(y))i>(y) = 0, (2.18) 

where 

^y) = 7^]2<^(y))' ( 2 - 19 ) 

J(y)^[r , (y)] 2 I(r(y)) + ^{r,y}, (2.20) 

The above procedure is called the Liouville transformation of \2. 1 6}) by the change of variables r = r(y) 
[Lioj . One can check using (I2.17[) , that the composition of two Liouville transformations, first by the 
change of variables r = r(y), then y = y(z), is the Liouville transformation by the change of variables 
r = r{y{z)). 



2.4 Natanzon's problem 

Let A(di, ■ ■ . , a.k', z, d z ) be a 2nd order differential operator depending on k parameters, such that its 
Bose invariant can be written as 

I(b 1 ,...,b n ;z) = b 1 I 1 (z) + ... + b n I n {z) (2.21) 

for some (at least one) n < k linearly independent functions I\, . . . , I n (independent of a±, . . . , a^), and 
some n numbers bi(a\, . . . , a^). Then one can solve the following problem (which we will call Natanzon's 
problem for A(z,d z )): 

Find all potentials V(r) such that the 1- dimensional stationary Schrodinger equation 

{-d 2 r + V(r) - E) <P{E, r) = (2.22) 

can be transformed to the equation given by A(a%, . . . , a^; z, d z ) for some a\(E), . . . , ak{E). We allow the 
following operations: 

(1) multiplication of both sides of the equation by some f(E,r); 

(2) substitution of (f)(E,r) := g(E,r)(j)(E,r) for some g(E,r); 

(3) change of coordinates r i— > z independent of E. 



This is the problem solved by Natanzon in the case of the hypergeometric equation and by Milson in the 
general case of hypergeometric type equations. In the following, we recall the construction they used, in 
a slightly generalized form. 

Let us consider an arbitrary equation of the form ()2.22|) . Obviously, its Bose invariant equals E—V(r). 
Clearly, the transformations (1) and (2) allow us to transform A(a\, . . . , a^; z, d z ) to its canonical form 

-d 2 z -I(bx,...,b n ;z). (2.23) 



G 



Thus (I2.22[) can be transformed to A{a-y, . . . , a/-; z, d z ) if 

I(h, . . . , b n ; z) = (r'(z)f (E - V (r(z)) + l{r, z}, (2.24) 

that is, if the two Bose invariants are related by a Liouville transformation. Using (|2.2ip we rewrite 
(|2~23|l as 

b x {E)h{z) + ... + b n {E)I n {z) = {r'{z)f (E-V(r{z)) + ^{r,z}. (2.25) 

By assumption, 6, depend on E and h(z) are independent of E. Thus, the dependence of 6, on E is 
linear. Therefore, by transforming linearly (b\, . . . , b n ) into (&i, . . . , 6 n _i, E), we can assume that 

I(bi, . . .,b n ;z) = 61/1(2) + . . . + 6 n _i/ n _i(z) + EI n (z), 

for some functions Ii(z), which are linear combinations of Ii(z). Then we can write 

EI n (z) + b 1 I 1 (z) + ... + b n _ 1 I n _ 1 (z) = (r'(z)) 2 (E-V(r(z)) + ^{r,z}. (2.26) 
Therefore, Natanzon's problem is solved by the following pair of equations 

V(h, 6 n _i,r(z)) := (r'(z)y 2 (bj^z) + ... + 6 B _ 1 7„_ 1 (z) - i{r, z} 

(r'(z)f = I n {z). (2.27) 
We can rewrite (|2.27l) in an equivalent way as 

V(pi, . . . , 6„_i,r) := (z'(r)) 2 (6 1 / 1 (z(r)) + ... + &„_ 1 /„_ 1 (z(r))) +i{z,r}, 

(z'(r)y 2 = I n (z(r)). (2.28) 

If desired, we can renormalize E by substracting a constant. 

Note that in practical applications the value of the above method depends on whether it is possible 
to invert the relation r <— > z expressing it in terms of standard functions. In the next sections, we 
will describe such solutions of the Natanzon problem for the three classes of equations described in the 
introduction. 

3 Schrodinger operators reducible to the hypergeometric equa- 
tion 

3.1 Hypergeometric equation 

The hypergeometric equation is given by the operator 

F(a, b; c; z, d z ) := z(l - z)d\ + (c-(a + b + l)z)d z - ab, (3.29) 

where a, b, c are arbitrary complex parameters. Its Bose invariant equals 

. , , (1 - a 2 - b 2 - 2ab)z 2 + 2(-2ab + ac + be - c)z + 2c- c 2 
2h(a,b;c;z) = 4 z 2 (l - z) 2 ' 
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It satisfies condition (|2.21[) for n = 3. The functions Ii(z), i — 1,2,3 can be taken to be 4z ^_ z ) 
4iH7-7P anc ^ 4z(i-z) • ^ * s na tural to demand that (r'(z)) 2 = I 3 (z) is a function of the form z 2p (l — z) 2q . 
It is easy to see that this gives the following possibilities for Iz(z): 

(1)7777^ (2)7-777— v (3) ' 



4z(l-*)' w 4z 2 (l-z)' w 4z(l-z) 2 ' 

< 4 > i^ibp (5) i' (6) i(iV (3 - 30) 

We will see that each of these ansatzes leads to exactly solvable potentials considered in the literature. 
Note also that the formula for the Schwarz derivative (j2.15[) will be handy. 

3.2 The Riemann equation 

It is well known that it is useful to consider the hypergeometric equation as a special case of the so-called 
Riemann equation. The Riemann equation is defined as the class of equations on the Riemann sphere C 
having 3 regular singular points. The following theorem summarizes the basic theory of these equations: 

Theorem 3.1 (1) Suppose that we are given a 2nd order differential equation on the Riemann sphere 
having 3 singular points Z\, z 2 , z 3 , all of them regular singular points with the following indices 

Zi : pi, pi, 

Z2 : P2,P2, 

Z3- p-i, Pi- 
Then the following condition is satisfied: 

Pi +Pl +P2 +P2 + P3 +P3 = 1- (3.31) 

Such an equation, normalized to have coefficient 1 at the 2nd derivative, is always equal to 

z,d z ct>{z) =0, (3.32) 




where 



p\ P \ pi pi z,d z I :=;> 2 - (^i-i + i^-i^ + ^-i)^ 



z — Z\ Z — Z<2 Z—Z3 



Pl P2 P3 



PlPl(zi-Z 2 )(zi-Z 3 ) P2p 2 (^2-Z3)(z2-Zl) , P3 P3 (^3 ~Zl ) (>3 ) 

' (z-z l y 2 (z-z 2 )(z-z 3 ) ' (z-z 2 )' 2 (z-z 3 )(z-z 1 ) ' (z-z 3 ) 2 (z-z 1 )(z-z 2 ) ' 

(2) Let z 1— > w = h(z) = ^frg. (Transformations of this form are called homographies or Mobius 
transformations). We can always assume that ad—bc= 1. Then 

h(zi) h(z 2 ) h(z 3 ) \ / zi z 2 z 3 

V I pi P2 p3 w,d w =(cz + d) 4 V\ pi p 2 p3 z,d z 



(3) 



Pl P2 P3 I \ Pl P2 P3 



(z- Zl )- X (z- z 2 ) x V pi p 2 p 3 z,d z \ (z~ Zl ) x (z- z 2 y 



Zl z 2 

V ' 



Zl 


Z2 Z 3 


Pl 


P2 P3 


Pl 


P2 P3 


Z3 


z,d z \ 


P3 




P3 





8 



Clearly, in all above formulas one of Zj can equal oo, with an obvious meaning of various expressions. 
For convenience we give the expression for the Riemann operator with 2:3 = 00: 

Z\ Z 2 OO 

Pi 92 P3 z,d z 
Pi h h 



' P1 + ~ P1 -1 pa+ft-r 



Z — Z\ Z — Z 2 

PlPl(zi-Z 2 ) p 2 p 2 (z 2 ~Zl) p 3 p3 



(3.33) 



(z - z 1 ) 2 (z - z 2 ) (z - z 2 ) 2 (z - zi) (z-z 1 )(z-z 2 ) 

The hypergeometric equation is a special case of the Riemann equation, since 

/ 1 00 

J^(a, b; c; z, d z ) — z(l — z)V \ a z,d z 

I 1-c c — a — b b 

= z(l - z)d z + (c - (a + b + l)z)d z - ab. 

Every Riemann equation can be brought to the form of the hypergeometric equation by applying (1) and 
(2) of Theorem O 

From Theorem 13.11 we also see that symmetries of the hypergeometric equation are better visible if 
we replace a, b, c with a, j3, p: 

a — c — 1 (3 = a + b — c, fi = a — 6; 
a= l±S±tJt } b= l±B±£±R : c = l + a. 

In fact, the new parameters coincide with the differences of the indices of the points 0, 1, 00: 

ot = pi-pi, P = p 2 - p2, P = Pa- pa- 
in the parameters a, /3, p, the Bose invariant of the hypergeometric equation has a more symmetric form: 

(l- a 2 )(l~2) + (l-/? 2 )2+( M 2 -l)(l-z)z 

2h{a,P,p,z) = 4z2(l - zf ■ 

We can summarize the relation between the hypergeometric equation and its canonical form by 
- z % _ z )%+h-ij: £±^±1; i + a] z, d z ) zS-i(l - z)-%~i 



- d * + ^ - ^ + 32 - l ) - { " 2 - 1} ^h) ■ (3 - 34) 



It will be natural to introduce the parameters 



\(a 2 -f3 2 ), 5:=\(a 2 +(3 2 ). (3.35) 



9 



In some cases the parameter k will be replaced by 

t := ^(a 2 -/3 2 ) =w. (3.36) 

We will describe two complex classes of exactly solvable potentials depending on two complex param- 
eters. Within each complex class there will be three real classes of exactly solvable potentials depending 
on two real parameters. 



3.3 Trigonometric Poschl-Teller potential 

In this subsection we consider Ansatz (1) of (|3.30|) . We set 



z = sin 2 — , which solves z' = z 2 (1 — z) 2 . (3.37) 



This leads to the Schrodinger equation 



where 



. l\ 1 COS T 

7 n 1" K 



4 



2 -2 

sin r sin r 



This potential was proposed and solved by G. Poschl and E. Teller |PTj . It is usually called the Poschl- 
Teller potential, sometimes also the Poschl-Teller potential of the first kind or the trigonometric Scarf 
potential. 

A natural real domain for this potential is ]0,7r[. If k,S are real, then the potential is real on this 
domain. 

Explicitly, the reduction of (|3 .38[) to the hypergeometric equation is derived as follows: 



-*t+i(i _ z )iH j- £±£_i£±i. i + a . Z}dz ^ z -f-i(i _ z y 

-z(l-z)V 



S. — L 

2 4 



1 oo 

a I 1 §_ i 1 a 

2 ' 4 2^4 2 

-9L J- I 4- I _E 

2 4 2^4 2 



-*(!-*)(«;+ (-L-^)a.) 

r \ 4^4sin 2 § V P 4^4cos 2 § 4 
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3.4 Hyperbolic Poschl-Teller potential 

We continue with (1) of (1530)) . We set 

2 = — sinh 2 ^, which solves z' = — (— z) 5 (l — z)^. (3-41) 
This leads to the Schrodinger equation 

^ 2 + < PT M + ^W) = 0, (3.42) 



where 



4 



coshr 



K- 



4 7 sinh 2 r sinh 2 r 

This potential was also proposed and solved by G. Poschl and E. Teller |PT] . In the literature it is known 
as hyperbolic, generalized Poschl-Teller potential, or the Poschl-Teller potential of the second kind. 

A natural real domain for this potential is ]0, oo[. If k,S are real, then the potential is real on this 
domain. 

To see that (|3.42[) can be solved in terms of the hypergeometric equation, we first repeat the compu- 
tations leading to (|3.40l) , and then set z = — sinh 2 § . 



3.5 Scarf potential 

We continue with (1) of (pT5UJl . We set 

I T T 11 

z = — — icosh — sinh — , which solves z = (— z)^ (1 — 2)2. (3.44) 
This leads to the Schrodinger equation 

(-a 2 + ^ s r (r) + ^^(r)=0, (3.45) 

where 

V s (r) * (s— 1 t sin ^ r 

V 4/ cosh 2 r cosh 2 r 

This potential was proposed and solved by F. Scarf [Sea] . In the literature it is often called the hyperbolic 
Scarf potential. 

A natural real domain for this potential is ] — oo, oof. If S, t are real, then the potential is real on this 
domain. 

To see that (|3.45p can be solved in terms of the hypergeometric equation, we first repeat the compu- 
tations leading to (|3.40l) , and then set z — | — i cosh | sinh | . 



3.6 More about Poschl-Teller and Scarf potentials 

Both kinds of Poschl-Teller potentials and the Scarf potential are real cases of the same complex case. 
To see this, consider eg. the hyperbolic Poschl-Teller potential as a function of the complex parameter 
r. It is holomorphic away from singularities at i7rn, n G Z. 
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For real 5, n, V$^ T (r) is real on iR and R + inn. On each halfline ]0, oo[+i7m and ] — oo, 0[+i7rn we 
obtain the hyperbolic Poschl- Teller potential. On each interval ]i7rn, i7r(n+l)[ we obtain the trigonometric 
Poschl- Teller potential. 

For real 5, in,, Vg% T (r) is real on R + in(n + |). On each of these lines we obtain the Scarf potential. 

Above, we used Ansatz (1) to derive Scarf and both kinds of Poschl- Teller potentials. Alternatively, 
one can use Ansatzes (2) or (3). To see this it is enough to consider Ansatz (3). In fact, we first repeat 
computations analogous to (|3.40j) : 



We set 



2 l\l-z B 2 ( 2 l\{l-z) 
a "4 — + "I 



z = tgh 2 -, which solves z' = z 2 (l — z), (3.46) 



obtaining the hyperbolic Poschl- Teller potential, 

z = -tg 2 ^, which solves z' = -(-z)^(l - z), (3.47) 
obtaining the trigonometric Poschl- Teller potential, or 

z = ctgh -, which solves z = z 2 (1 — z), (3.48) 

obtaining the Scarf potential. 

3.7 Manning-Rosen potential 

Let us consider Ansatz (3) of (|3.30p . We set 

z = - — 2r , which solves z' — 2z(z — 1). (3.49) 
This leads to the Schrodinger equation 

("# + O r )+ 5 )^)=0, (3.50) 

where 

V MR (r) ■= _ _ 1^ 1 

^coshr U 4y G osh 2 r' 

This potential was proposed and solved by M. F. Manning and N. Rosen [MR] . In the literature it is also 
called the Woods- Saxon potential |SWj (for instance in [Fliij ) . also the hyperbolic Rosen-Morse potential. 

A natural real domain for this potential is ] — oo, oof. The potential is real if r £] — oo, oof and k, j3 2 
are real. 
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Here is a derivation of (I3.50[) from the hypergeometric equation: 

-4z 1+ t(l - z)i+fjT( g±g±H±l , £±£_i£±i . i + a; Zi d,)z-%(1 - z)- 



1 oo 

2 2 ^ 2 



-4z 2 (l-z) 2 ^<9 2 + - Y^i) 3 *) +a 2 {l-z) + p 2 z-{fi 2 ~\)z(l-z) (3.51) 



1 



3 - + 5 + ^1^27 - 1 ) (l +e 2r)2 - ( 3 - 52 ) 



3.8 Eckart potential 



We still consider Ansatz (3) of (|3.30p . 

We set 

z = -5-, which solves z' = 2z(l — z), (3.53) 

1 — e lr 

and use (|3.35p . This leads to the Schrodinger equation 

K 2 + O r ) + 5 )^) = °. ( 3 - 54 ) 

where 

E _ coshr / [i 2 l\ 1 
K ^ (2/) := ~ K ihm7 + IT ~ ij TnnV 

This potential was proposed and solved by C. Eckart [Eckj . In the literature (for instance [Flii ) it is 
also called the Hulthen potential [Hulj . sometimes also the generalized Morse potential, because of its 
similarity to the Morse potential, see Subsect. 14.41 

A natural real domain for this potential is ]0,oo[. If K, (3 2 are real, then the potential is real on this 
domain. 

To derive the Eckart potential, we first repeat the computations leading to (|3.51[) . and then set 

2 = l-e- 2 ^" 

3.9 Rosen-Morse potential 

Once again, we consider Ansatz (3) of (|3.30p . 
We set 

z= l ~2ir' which solves z' = 2iz(l — z). (3.55) 

This leads to the Schrodinger equation 

(-9 r 2 + K™(r) - 6) J>(r) = 0, (3.56) 

where 



sin r \ 



• 2 ' 

sm r 



cosr ( ir 1\ 1 
~~ 4, 

This potential is known as the Rosen-Morse potential, also the trigonometric Rosen-Morse potential 
(altough this name is widely used in the literature, we were unable to explain decisively its origin). 
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A natural real domain for this potential is ]0,7r[. If r, are real, then the potential is real on this 
domain. 

To derive the Rosen-Morse potential, we first repeat the computations leading to (I3.51[) . then set 
z - T=^- 

3.10 More about Manning- Rosen, Eckart and Rosen-Morse potentials 

The Manning-Rosen, Eckart and Rosen-Morse potentials are all real cases of the same complex case. To 
see this, consider eg. the Eckart potential as a function of the complex parameter r. It is holomorphic 
away from singularities at i7rn, n £ Z. 

For real k, [J. 2 , V^Ar) is real on the lines R + On each line R + i(m + |)7T, m € Z, we obtain 
the Manning-Rosen potential. On each halfline ]0, oo[+i7rm and ] — oo, 0[+i7rm, we obtain the Eckart 
potential. 

For real in,^ 2 , V^F„(r) is real on iR. On each interval ]i7rm, i7r(m + 1)[ we obtain the Rosen-Morse 
potential. 

Ansatzes (5) and (6) lead to the same classes of exactly solvable potentials as Ansatz (4). To see this 
it is enough to consider Ansatz (5). In fact, we first repeat computations analogous to (|3.51|) : 

-48-^(1 - z) 2 +l,F( g±g±g+i , i + a . z% d z )z-^{l - z)S 




-A{l- Z yid 2 z - — 



- I ^ 1 -n. 

z 



+ (a 2 - 1)— 5- + /3 2 - - M • 

z z z z 



We set 



z = 1 - e 2r , which solves z = -2(1 - z), (3.57) 
obtaining the Eckart potential, 

z = l + c 2r , which solves z' = 2(1 - z), (3.58) 
obtaining the Manning-Rosen potential, 

z = l-e i2r , which solves z' = -2i(l - z), (3.59) 
obtaining the Rosen-Morse potential, 

4 Schrodinger operators reducible to the rescaled confluent equa- 
tion 

One of basic exactly solvable equations is the confluent equation, given by the operator 

T{a; c; z, d z ) := zd 2 + (c - z)d z - a. (4.60) 
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It is convenient to consider (|4.60l) in parallel with the equation given by the operator 

F(a, b;-;z, d z ) := z 2 <9 2 + (-1 + (1 + a + b)z)d z + ab. (4.61) 

The equation given by (|4.61[) is sometimes called the 2-F0 equation. Note that 

z a J r (a, b; — ; z, d z )z~ a — wJ-(a; 1 + a — 6; w, d w ), w — —z^ 1 ^ z = —w . 

Hence the 2-F0 equation is equivalent to the confluent equation. The relationship between the parameters 
is 

c = 1 + a — 6, & = 1 + a — c. 

Another exactly solvable equation that we will consider in this section is sometimes called the qFi 
equation. It is given by 

F(c;z,d z ) := zd 2 + (c - z)d z - a, (4.62) 
and is equivalent to the Bessel equation. 

Clearly, the confluent, 2-F0 and qF\ equations belong to the class of hypergeometric type equations. 
The basic equation of this section will be (jl.5[) , given by 

T{a; c; 7; z; d z ) := zd z + (c — jz)d z — a. (4.63) 

It will be called the reseated confluent equation. Note that in the case 7=1 (|4. 63[) coincides with the 
confluent equation. If 7 7^ 0, (|4.63[) can be reduced to the confluent equation by rescaling (and hence 
also to the 2-F0 equation). If 7 = 0, (|4.63[) coincides with the 0-F1 equation. 
The Bosc invariant of the rescaled confluent equation equals 

T , , N -7 2 z 2 + 2(c 7 - 2a)z + 2c - c 2 

ih(a,b;r,z) = 7-5 ■ 

4z z 

4.1 Symmetries of the rescaled confluent equation 

Let us first describe symmetries of the rescaled confluent equation. 

z c - 1 J r (a- 1 c;j- 1 z,d z )z 1 - c = F(a - cry + 7; 2 - c; 7; z, d z ); (4.64) 
e~ lz J 7 (a; c; 7; z, d z )e JZ = — T{c — a; c; 7; w, d w ), z = —w; (4.65) 
Besides, the scaling acts as follows: 

T(a; c; 7; z, d z ) = 7^(0/7; c; 1; w, d w ), w = -yz. 
It is convenient to introduce new parameters a, v. 

a = ' 2 ' — , c = 1 + a: 

a = c — 1 = a — 6, v — cy — 2a = I — a — 6. 
In the new parameters the Bose invariant of the rescaled confluent equation has a more symmetric form: 

7 2 v 1 — a 2 

iMw;*) = -~ + Yz + ^z^- 

Thus the starting point for the further analysis will be the equation 

-«4-5 + (T-i) ?)«"-■ <« 6 > 

We will describe 3 classes of Schrodinger operators solved using the confluent equation corresponding 
to three obvious choices for (r'(z)j : 

« I> ( 2 ) (3) (4-67) 
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4.2 Hydrogen atom 

We consider (1) of (|¥T571i . We set 

z = 2r. (4.68) 

This leads to the Schrodinger equation 

{-d 2 r + V a , v (r) + 7 2 ) 4>(r) = 0, (4.69) 

with the potential 

T , v fa 2 1\ 1 

: =-7 + U-4j^ 

(|4.69p is the radial part of the Schrodinger equation for the hydrogen atom. 

]0, oo[ is a natural real domain for this equation. If v, a 2 are real, then so is V atV on ]0, oo[. 
The derivation of (|4.69[) from (|4.66[) is immediate: 

z 4 22 U 4/z 2 

^2 2^ f°? 1\ 1 

4.3 Rotationally symmetric harmonic oscillator 

We consider (2) of (j4~67|) . We set 

z = r 2 , which solves z' = 2^/z. (4.70) 
This leads to the Schrodinger equation 

{-d 2 r +V a , 1 {r)-2v)4>{r) = 0, (4.71) 

with the potential 

V a „{r) :=~, 2 r 2 + 



- ' <a 2 --^ ' 



4y r 2 ' 

]0,oo[ is a natural real domain. For real j 2 ,a 2 , the potential is real on ]0,oo[. By scaling, the only 
different real cases are 7 2 = 1,0,-1. 

For 7 2 = 1 the equation is the radial part of the rotationally symmetric harmonic oscillator. For 
7 2 = the equation is the radial part of the Helmholtz equation. 

Here is an explicit derivation of (|4.T1[) from (|4.66l) : 



Az^(-d 2 + 1 --- 1 

\ z + 4 2z ' V 4 47 z 2 

= -4zd 2 - 28 z + 7 2 z - 2v + [a 2 - ~) i 

= -a r 2 + 7 V-2,+ (^-i)i. 



z< 
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4.4 Morse potential 

We consider (3) of We set 

z — e~ r , which solves z' — —z. (4.72) 
This leads to the Schrodinger equation 

v 2 
4 



^ + K, 7 (r) + ^U(r)=0, (4.73) 



where 

V v , 1 {r) = ^~ 2r -^- r - 

] — oo, oo [ is a natural real domain. For real v, 7 2 , the potential is real on ] — oo, oo[. By translation, 
the only different real cases are j 2 = 1,0, — 1. 

Here is an explicit derivation of (|4.73l) from (|4.66l) : 

^(-* + t-s + (t;i)f)^ 

= -z 2 dl-zd z + ^- V - Z+ ^ 
r + 4 2 4 ' 

5 Schrodinger equations reducible to the translated harmonic 
oscillator 

The last family of exactly solvable Schrodinger equations that we describe is (|1.6|) . given by 

- d 2 + 6» 2 z 2 + + A. (5.74) 

We will call (|5.74[) the translated harmonic oscillator equation. It is already in the canonical form. 

If 9 2 ^ 0, then (|5.74j) is just the translation of the Schrodinger equation for usual harmonic oscillator. 
It is then equivalent to the Hermite equation, given by 

G(a,y,d y ) :=(d 2 -2yd y ~2a). 

In fact, we have 

2 2 

9e~^G(a,y,d y )e^ 
= e(-d 2 y +y 2 + 2a-l) 

= -d 2 z +e 2 z 2 + pz + ^ + 2a-l, 

where y = \f§ (z + . 

For = 0, p 7^ 0, (|5.74|) is equivalent to the Airy equation 

(a 2 + y)i>(y) = 0, 
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which in turn is equivalent to a special case of the o-Fi equation: 

-d 2 + pz + \ 
= pH-d 2 y +y) 

= (p3)3w3j- Q; Wi d w , . 

where we set y — p's z + p~~s A, and w = 3~ 2 y 3 . 

We will describe 3 classes of Schrodinger operators that can be reduced to (|5.74[) corresponding to 



three obvious choices for 



(1) 1, (2) z, (3) z 2 . (5.75) 



5.1 Translated harmonic oscillator 

We consider Ansatz (1) of (|5.75[) . which corresponds to the most obvious choice of the energy, that is —A. 
We rename the variable z — r. This leads to the Schrodinger equation 

(-d 2 + V 0tP (r) + A) <j>{r) = 0, (5.76) 

where 

Vgjr) =6 2 r 2 +pr. 

5.2 Special potential I 

We can choose the energy to be —p, that is Ansatz (2) of f|5 . T5|) . This corresponds to the substitution 



z 



3r\ 3 

— ) . which solves z — z _3 , 
2 



and leads to the Schrodinger equation 

{-d 2 + Vg.x(r)+ p)<t>(r) = 0, (5.77) 



where 



In fact, 



2 ' 



z^- 1 (-a% + e 2 z 2 + pz + x) 

Note the following intriguing feature of the above potential: the coefficient at r~ 2 is fixed and one 
cannot change it by rescaling the variable r. 
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5.3 Special potential II 

We can choose the energy to be — 9 2 , that is Ansatz (3) of (|5.75|) . This corresponds to the substitution 

z = (2r) 5 , which solves z' = z -1 , 
and leads to the Schrodinger equation 

(-d 2 + V pA (r)+8 2 )cj)(r) = 0, (5.78) 

where 

I \ P ^ 3 1 

V pM r ) = TT-T + 7T ~ TZ~- 

(2r)s 2r 16 r z 

In fact, 

z 5-2 (_^2 + ^2 + pz + \) J,— * 

Again, the coefficient at r -2 is fixed and one cannot change it by rescaling the variable r. 
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